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Using post-Newtonian equations of motion for fluid bodies that include radiation-reaction terms 
at 2.5 and 3.5 post-Newtonian (PN) order (0[(i>/c) 5 ] and 0[(v/c) 7 ] beyond Newtonian order), we 
derive the equations of motion for binary systems with spinning bodies, including spin-spin effects. 
In particular we determine the effects of radiation-reaction coupled to spin-spin effects on the two- 
body equations of motion, and on the evolution of the spins. We find that radiation damping causes 
a 3.5PN order, spin-spin induced precession of the individual spins. This contrasts with the case of 
spin- orbit coupling, where we earlier found no effect on the spins at 3.5PN order. Employing the 
■ equations of motion and of spin precession, we verify that the loss of total energy and total angular 

' momentum induced by spin-spin effects precisely balances the radiative flux of those quantities 

calculated by Kidder et al. 
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INTRODUCTION AND SUMMARY 



The relativistic effects of spin may play an important role in the inspiral of compact binary systems, particularly 
involving black holes, and may have observable effects on the gravitational-wave signal emitted. Spin-orbit and 
spin-spin coupling leads to precessions of the spins of the bodies and of the orbital plane, the latter effect resulting 
in modulations of the amplitude of the gravitational waveform received at a detector. Furthermore, spin effects 
contribute directly to the gravitational waveform, and to the overall emission of energy and angular momentum from 

■ the system. 

f-*. | In the post-Newtonian (PN) approximation to general relativity, the effects of spin have been derived at various 

■ levels of the PN approximation. Formally, spin-orbit and spin-spin couplings begin to affect the equations of motion at 
""^ | the first post-Newtonian order, since they behave as S-L/mr 4 ~ Si -S2/mr 4 ~ (mRV)(mrv) / mr A ~ (mRV) 2 jrar^ ~ 

Q" 1 (m/r 2 )e, where m, v, r and L denote mass, orbital velocity, separation and orbital anglar momentum, respectively, R 
and V denote the body's size and rotational velocity, S denotes spin or rotational angular momentum, and e ~ v 2 ~ 
(3Tj, V 2 ~ m/r ~ mjR denotes the standard small "bookeeping parameter" of post-Newtonian theory (we use units in 
which G = c = 1) 1]. Spin evolution effects can also be seen to be lPN-order effects. Indeed, the 1PN effects of spin 
. ! have been derived by numerous authors from a variety of points of view, ranging from formal developments of the 
GR equations of motion in multipole expansions 0, Q , to post-Newtonian calculations [H , to treatments of linearized 
GR as a spin-two quantum theory [H, Q. For a review of these various approaches, see Q- The effects of spin on 
the gravitational waveform and on the energy and angular momentum flux were worked out by Kidder et al. [1, 0] . 
Post-Newtonian corrections of the leading spin terms have also been analysed [l(| EH ■ 

In Paper III of this series (T3 |. we derived, from first principles, the leading effects of spin-orbit coupling in the 
equations of motion at radiation reaction order, specifically at 3.5PN order, or 0(e 7 / 2 ) beyond Newtonian gravity. We 
also showed explicitly that radiation reaction had no effect, via spin-orbit coupling, on the individual spins themselves. 
In this paper, we extend this analysis to spin-spin coupling. As before, the leading contributions occur at 3.5PN order. 

We use the hydrodynamic equations of motion derived through 3.5PN order in Papers I and II [TH, HU, and 
calculate the equations of motion and spin evolution for two spinning, finite-sized bodies. We restrict our attention 
to contributions that involve the products of the two spins. To this end, for each body A, we decompose velocities 
into a center-of-mass part and an internal (rotational) part according to v = + v, and expand all gravitational 
potentials about the center of mass of each body using a similar decomposition, x = + x, and retain only terms 
that contain the product (mxv)i(mxv)2 ~ S%S2- We do not keep terms proportional to the squares of individual 
spins; these represent another class of spin effects that will be studied elsewhere. 

Adopting a specific definition of "proper spin" Sa, as defined in Paper III [see Eq. (|2.15p ], we find the two-body 
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equations of motion 



a = — + apN + apN-so + apN-ss 



&2.5PN + a3.5PN + a3.5PN-SO + a3.5PN-SS 



where a = a x — a 2 is the relative acceleration. The 1PN spin-spin terms are standard, and are given by 

3 

apN-ss = 4 [n(5i • S 2 ) + Si(n • S 2 ) + S 2 (n • Si) - 5n(n • Si)(n • S 2 )] , 

and the 3.5PN spin-spin contributions, derived in this paper, are given by 

1 f \( 7 , 541 m\ / , 9 1961 m\ , 

287f 2 - 99v 2 + — )r(Si ■ S 2 ) - 2646r 2 - 7Uv 2 + r (n • Si) (n • S 2 ) 

or/ \ or) 



(1.1) 



(1.2) 



a3.5PN-SS — -g < ii 



1029f 2 - 123?; 2 
171 



629 m 

To"7 



(n • Si) (v • S 2 ) + (n • 5a) (v -Si))- 336f (v • Si) (v • S 2 ) 

m\ / 1038 m 

-^v 2 - 195f 2 - 67— (Si • S 2 ) - 174w 2 - 1386f 2 I (n • Si) (n • S 2 ) 

5 r / V or 



438f (n • 5i) (v • S 2 ) + (n • S 2 ) (v • Si) + 96 (v • Si) (v • S 2 ) 



27 2 75. 2 509m , 
10 V ~Y T ^^7f7 )( ( v ' 5 2)5i + (v-Si)S 2 



15 



77 



199 m 



f^ + f^ + lif ) f ( (n-5 2 )Si + (n-Si)S 2 



(1.3) 



where x = xi — x 2 , r= |x|, n = x/r, m = mi+ m 2 , fi = mim 2 /m, v = vi — v 2 , r = dr/dt = n • v, and Ln = xxv. 
The PN, PN spin-orbit, and 2.5PN contributions are standard (see Paper III Eqs. (1.2) and (1.5) for the formulae), 
the 3.5PN terms were derived in Paper II, Eq. (1.3d), and the 3.5PN spin-orbit terms were derived in Paper III, Eq. 
(1.6). 

The equations of spin evolution are given by 



Si — (Si)pn-SO + (5l)pN-SS + (5l)3.5PN-SS : 



(1.4) 



where the PN spin-orbit terms are standard (see Paper III, Eq. (1.3)). The PN spin-spin terms are also standard, 
given by 



1 



lJPN-SS 



S 2 - 3 (n • S 2 ) n x Si 



There is no 3.5PN spin-orbit contribution (Paper III), but we find a 3.5PN spin-spin contribution, given by 

jim ( 2 



(Si) 



3.5PN-SS 



(v-S 2 ) + 30f(n-S 2 ) (n x Si) , 



(1.5) 



(1.6) 



with the equations for S 2 obtained by interchanging the spins. 

As a check of these results, we verify explicitly that the loss of total energy and total angular momentum (including 
both orbital and spin) implied by these equations matches the energy and angular momentum radiated in gravitational 
waves, as calculated by Kidder et al.\^, Q. 

In Paper III, we found that spin-orbit contributions to radiation reaction had no effect on the proper spin of each 
body, i.e. (Si)3.5pn-so = 0> an d we argued that this made sense, given that a spinning, axisymmetric body should not 
couple to gravitational radiation. Here, however, when the coupling between the two spins is taken into account, there 
is a non-trivial radiation-reaction effect on the spins. Nevertheless, the effect is a pure precession; the magnitude 
of the spins is unaffected. Furthermore, if either of the spins is aligned with the orbital angular momentum (i.e. 
perpendicular to v and n), the other spin is not affected by radiation reaction. 

These equations of motion do not impose any limitations on the orbits. In particular, they can be used to evolve 
the quasi-circular inspiral orbits that are typical of those considered as sources of gravitational radiation detectable by 
ground-based later-interferometric detectors of the LIGO- VIRGO type, as well as highly eccentric orbits of extreme 
mass ratio systems that are relevant for the proposed space-based detector, LISA. 

The remainder of the paper presents details. In Section|TT]we derive the equations of motion to PN order, including 
spin-spin terms, and show that no spin-spin effects occur at leading radiation-reaction, or 2.5PN order. This section 
illustrates some basic features of the technique of obtaining the spin effects from the hydrodynamical equations. In 
Section IlIII we move to 3.5PN order, where the spin-spin radiation reaction effects first appear. Section llVI presents 
concluding remarks. 
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II. POST-NEWTONIAN AND 2.5PN EQUATIONS OF MOTION AND SPIN EVOLUTION 

A. Foundations 

As in Paper III 1^], we analyse a binary system consisting of balls of perfect fluid that are sufficiently small 
compared to their separation that tidal interactions (and their relativistic generalizations) can be ignored, but that 
are sufficiently extended that they can support a finite rotational angular momentum, or spin. At Newtonian order, 
the result is essentially trivial: the equation of motion for body 1 is d 2 Xi/dt 2 = —m^x/r 3 + 0(mR 2 /r 4 ), where R 
is the characteristic size of the bodies. Spin plays no role whatsoever, because the Newtonian interaction does not 
depend on velocity. But at post-Newtonian order, there are velocity-dependent accelerations of the schematic form 
mv 2 /r 2 , and thus, taking into account the finite size of the body and expanding about its center of mass, we expect 
to find acceleration terms of the form (mVR)(mVR)' /mr A ~ SiS^/mr 4 . However, the combination of finite size and 
spin introduces an ambiguity in the definition of the center of mass of each body. This has given rise to the concept of 
"spin supplementary condition" (SSC), a statement about which center of mass definition is being used; this concept 
is discussed in Paper III, Appendix A. It turns out that this is an issue only for spin-orbit effects; the choice of SSC 
or of center of mass has no effect on spin-spin effects at PN or at 3.5PN order. 

We will define centers of mass and spins provisionally using the "conserved" , or baryonic density, given by 

p* = p^u\ (2.1) 

where p is the mass energy density as measured by an observer in a local inertial frame momentarily at rest with 
respect to the fluid, g is the determinant of the metric, and u° is the time component of the fluid four-velocity. 
Assuming that p is proportional to the baryon number density, then conservation of baryon number leads to the 
useful exact continuity equation 

dp*/dt + V • (p*v) = , (2.2) 

where v l = u % /u is the ordinary (coordinate) velocity of the fluid. (Greek indices range over spacetime values 
0, 1, 2, 3, while Latin indices range over spatial values 1, 2, 3. Henceforth, spatial vectorial quantities will be handled 
using a Cartesian metric.) The baryonic mass, center of mass and baryonic spin of each body in our system are 
defined to be 



m A = I p*d 3 x, (2.3a) 

A 



xa = m 



1 / p*xd 3 x, (2.3b) 

J A 



S A ee / p'xxvft, (2.3c) 

J A 

where x = x — xa and v = v — . We also define a two- index spin quantity 

S% ee 2 / p*x [l v l] d 3 x, 



A 



e ijk S k A , S\ = \<# h S* , (2.4) 



where [] around indices denotes antisymmetrization. With these definitions, the baryonic mass rriA is constant, and 
the velocity, acceleration and rate of change of spin of body A are given by 

va = m^ 1 / p*vd 3 x, (2.5a) 

J A 

sla = m ~A / P*&d 3 x, (2.5b) 

J A 

dS A /dt = / p*xxad 3 x. (2.5c) 
J A 

Notice that, by virtue of the definition of center of mass, the "bar" can be dropped from the acceleration in Eq. 
(|23c| . 
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B. Baryonic equations of motion and spin evolution 



We begin by working to 1PN and 2.5PN order, reproducing the standard 1PN formulae for spin-spin interactions, 
and establishing some results that will be useful when we go on to 3.5PN order. Since we are only interested in 
radiation-reaction aspects of spin, we can ignore the 2PN terms in the equations of motion; these produce only 
conservative PN corrections to the spin equations of motion 

We use the hydrodynamic equations of motion derived in Paper II, Eqs. (2.23), (2.24a) and (2.24c), with all 
quantities expressed in terms of the conserved density p* . They are given by 



d x % jdt = U' % + a s PN + fl2.5PN i 
where U is the Newtonian potential, and where 



(2.6) 



a 



PN 



v^U' 1 - IrVT •' - 4UU' 1 - iv l U + W l + 8v j V [i ' j] + - $ 9 l + -X< 



(I 



o (5) (5) (4) (3) 4 (3) (3) 

_ pN = -stUp* --S ij l kk ) + 2v? T j +2U- j /'•' +-U' i l kk ~X^ k l jk , 



(2.7a) 
(2.7b) 



where commas denote partial derivatives, overdots denote partial time derivatives, and (n) above quantities denotes 
the number of total time derivatives. The potentials used here and elsewhere in the paper are given by the general 
definitions 



X(f) = J p*(t,x')f(t,x')\x-x'\d 3 x 



(2.8) 



with specific potentials given by 



U = E(l), V* = £(«*), 

$i = iV), $ 2 = s(c/), 

X = X(l), X 1 =X{v l ). 



(2.9) 



The multipole moment of the system I y , as well as additional moments, J 1 ^ , J l i , and A4 t3 , that will be relevant 
at 3.5PN order, are defined in Paper III, Appendix C (see also Eq. (|3.2[) ); note that there are no explicit spin-spin 
terms in X y , to the PN order considered. 

We now multiply the equation of motion (|2.6[) by p* and integrate over body 1, expressing the variables x and v as 
x = x_4 + x and v = + v, where A = 1, 2, depending on the body in which the point lies. To get the acceleration 
of body 1, we divide the result by mi. We use Eqs. (|2.3|) and (|2.5p to simplify where possible. We expand the various 
potentials in powers of x/r, and keep only terms proportional to the product of v x x for one body with iixi for the 
other body. 

In Paper III, we also kept internal terms proportional to v 2 and used virial relations derived in Paper III, Appendix 
E to simplify expressions dependent on the internal structure of each body. While the use of those virial theorems 
generated spin-orbit terms at PN order, it turns out that they generate no spin-spin terms at this order. We will deal 
with the effect of virial theorems on 3.5PN spin-spin terms in Sec. IIIII 

The Newtonian term gives = — m 2 x l 12 /r 3 , where, in this paragraph, to avoid confusion, we denote 
r = |xi2 1 , and n = X12/V. The only PN terms in Eq. (|2.7a[) that can have a v in one body and a v in the other 
body and that could therefore lead to a spin-spin effect are the terms — 3v l U and Sv^V^^ . Keeping only the relevant 
terms, the first of these gives, for example, 



3 — / p*v i Ud 3 x 
mi J x 



-3— / p*(v\ +v z )d 3 x 
mi J x 



x J p*'{J 2 +vP) 



% + (x-x') k V k ( X \ 



+ l^-x') k {x-x') l V k V l 



12 



Anii 



s kl s l 2 J v k v l pf , 



d 3 x' 



(2.10) 
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where unprimed (primed) barred variables are in body 1(2). We also assume that each body is in stationary equilib- 
rium, with if = (d/dt) J A p*x i x j d 3 x = 0, so that J A p*x i vH s x = S l j/2. 

In the combination of PN terms 4V l + \X l in Eq. (|2.7a|l . the time derivatives generate accelerations inside the 
potentials. To the order needed for our purposes, we must therefore substitute the Newtonian and 2.5PN continuum 
terms for those accelerations and carry out the same procedures for the integrals as described above. However, while 
this produces spin-orbit effects (Paper III), it produces no spin-spin terms. 

The resulting 1PN spin-spin contribution to the equation of motion for body 1 is 

(ai)pN-ss = — —j [sVS%n k + S^SfV + sFsgn* - 5SpS l jn i n k n l ) . (2.11) 
mir 4 V / 

Substituting Eq. (|2.4[) and calculating the relative acceleration gives 

apN-ss = - A t n ( S i ' S 2) + Si(n • S 2 ) + S 2 (n • Si) - 5n(n • Si)(n • S 2 )] , (2.12) 

where the spins here are baryonic spins. 

Turning to the 2.5PN terms, Eq. (|2.7bp . the integrations lead to no explicit spin-spin terms, so that the 2.5PN 
relative acceleration terms are given by 

o (5) i (5) (4) 1 (3) (3) 

(a% 5PN = -x 3 (l 11 --S 11 l kk ) + 2v* T 1 ~-^n l l kk -3^n l n 3 n k V k . (2.13) 
5 6 A r A r £ 

However, when we work at 3.5PN order, even though the multipole moments themselves contain no explicit spin- 
spin terms, time derivatives acting on them will produce spin-spin contributions via the PN spin-spin terms in the 
equations of motion. 

We calculate the evolution of the spin in a similar manner. Starting with dS\/dt — e 13 L p*x 1 a k d 3 x, we expand 
about the baryonic centers of mass, keeping only terms that depend on a product of x, x v for each body. At 1PN 
order, the only term in Eq. (|2.7a[) that contributes is Sv^V^. The result, at 1PN order is 

(Si)pn-ss = fs a - 3(n ■ S 2 )n) x Si , (2.14) 
where again these are baryonic spins. At 2.5PN order, there is no spin-spin contribution. 



C. The proper spin 



In Paper III, we defined the proper spin of each body to be 

/ i \ 1 (3) (3) 

S[ = sUl + -v\ + - -[vi x ( Vl x Si)]* - S{ 1" +S{ V 1 , 

S 2 = (1^2). (2.15) 

The post-Newtonian corrections in Eq. I|2.15p arise from transforming the baryonic spin from our coordinate frame 
to a suitable inertial frame comoving with the spinning body. The 2.5PN terms involving time derivatives of X 13 arise 
from the fact that the equations of motion at 2.5PN order may be written in various gauges, such as Burke-Thorne 
gauge [l5j (in which the 2.5PN radiation reaction terms in the acceleration are given by a\ 5PN = |a;-'d 5 2" <lJ> /(it 5 ), or 
Damour-Deruelle gauge [S[l7| (the gauge used in this paper). Including the 2.5PN terms as in (|2.15|) is equivalent 
to defining our spins in the Burke-Thorne gauge. In any case, such quantities as angular momentum and energy are 
well-defined only up to the order at which they conserved, and one is free to add 2.5PN and 3.5PN order terms to 
them without affecting their fundamental conserved properties; including the 2.5PN terms in (|2.15j) has the property 
that radiation reaction effects in the proper spin do not appear (if they appear at all) until 3.5PN order. With 
this definition, the proper spins 5,4 also satisfy the standard spin-orbit precession equations, Paper III, Eq. (1.3). 
Transforming from our baryonic spin to the proper spin will generate some spin-spin terms at 3.5PN order. 
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D. Conserved total energy and angular momentum 

The Newtonian and PN spin-spin terms in the equation of motion (|2. 12|) . and the PN spin-spin terms in the spin 
precession equation (|2.14p . together imply conservation of the total energy and angular momentum of the system, 
given to Newtonian and PN spin-spin order by 

/ 1 m\ 1 

E = p ( - v 2 - — j - — [(5i ■ 5 a ) - 3 (n • Si) (n • 5 2 )] , (2.16a) 

J = pL N +S , (2.16b) 

where 5 = 5i + 52- Notice that there is no PN spin-spin contribution to the total angular momentum. In Eq. 
(|2.16[) we have converted the baryonic spins to the proper spins; the PN corrections in Eq. (|2.15[) do not introduce 
new spin-spin contributions to E or J to the required order, and the 2.5PN corrections can always be dropped as 
meaningless terms that have no effect on the conserved quantities. 

In Sec. IIIICI we will use these expressions together with the 3.5PN equations of motion to compare E and J with 
the corresponding fluxes of radiation to infinity. 



III. 3.5PN EQUATIONS OF MOTION AND SPIN EVOLUTION 
A. Equation of motion 

To obtain the 3.5PN contributions to the equations of motion including spin-spin coupling terms, we take the 
3.5PN fluid expressions shown in Eq. (D4) of Paper III, multiply by p* , and integrate over body 1. We follow the 
same procedure as in Sec. Ill B| expanding potentials about the baryonic centers of mass of the bodies, keeping only 
spin-spin terms (terms involving products of x x v for one body with that for the other body) . Most terms in Eq. (D4) 
of Paper III make only point-mass or spin-orbit contributions; the only terms that can possibly produce non-trivial 
spin-spin terms in either the acceleration or the spin evolution are: 

(3) 1fi (3) (3) 1fi (5) 

Sal , PN = -8v k V k ' j T ij +— v j V^ I kk +8(v j V k ' i - v l X^ jk ) l jk x j v k e qjk J qi 

2 (5) 
+ — (2(v • x)e qlk - 2x l v ] e qik + hx 3 v l t q3k + I2x j v k e qij + ix k v j e qij ) J qk 

2 (4) 2 (3) 

+-(4v j v k e qij - v 2 e qik ) J qk -~(^t/' j + 2v j U' i + + X'^)e qjk J qk 
9 9 

1 O) (5) (5) i (4) „ (4) 

+ _ u J( e ?j'fc jqik _ e qik jqjk _ £TO jqkk\ _ _ y i M kk n + _ v j M ijkk ^ (3 ^ 

15 6 3 



where, to the required order, 

^12^12 ! 



1* ^ — pX-^ 2 X^- 



J** = -vSmLiixiz - ^(3AV 12 - S* A ■ x 12 ) , 
J» k = r,m(l - 3 V )V N x{ 2 x k 12 + ^x{ 2 x\ 2 - £ ■ x u ^) , 
M ijkk = vm{1 _ 3v y ^| 2 ^ 2 _ ^.n\ 2 n{ 2 ) - ± V m 2 r(n\ 2 n{ 2 - 3^) - 2r?(x 12 x OHl > (3-2) 

where 5m = mi-m2, r\ = p/m, £ l = (m2/m{)S\ + (rni/m2)S 2 , and A 1 — m(S 2 /m2~ S\/m\). Spin-spin contributions 
come from terms such as v^V k ' 1 d 3 I^ / dt 3 , with a v in body 1 and a v in body 2 together with suitable s's. They also 
come from terms involving the current moments J qk , where a single spin generated by the prefactor (eg. asVe 9 - 7 , or 
v z U " ,J 'e* 7 ) is multiplied by the spin of the other body that appears in J qk . The terms involving d 4 J qk /dt 4 d 5 J~ qpk / dt 5 
and d 4 Ai pqkk / dt 4 don't generate spin-spin terms in the equation of motion (no free x to go with a velocity), but do 
generate terms in the spin evolution. 

In addition, when the prefactor of d 3 J qk /dt 3 is integrated over body 1, it yields a prefactor given by — (2/9)(47i^"'- ) — 
3/C^), where 

^ s i A i A p * p ^\^? Mx '' s l i A p * p * - (x " x ix - x y {x ~ — Mx> ■ (3 - 3) 
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However, a virial theorem derived from the requirement that = gives, to the required PN order, 




(See Paper III, Appendix E, for a discussion of virial relations.) This spin term, combined with those generated 
directly by the potentials, and multiplied by the appropriate spin term in d 3 J qk / 'dt 3 , gives a 3.5PN spin-spin term. 
This is the only place where the virial theorems play a role in spin-spin effects. 

In addition, the combination of 1PN terms 4V l + ^X' % in Eq. (|2. 7a)l . will generate accelerations whose 2.5PN 
terms will produce 3.5PN contributions; however these produce no spin-spin terms. We must also re-express the 1PN 
spin- spin terms of Eq. p. lip in terms of the proper spin of Eq. (|2.15p : the 2.5PN contributions there will generate 
3.5PN spin-spin terms in the equation of motion. Finally, in the 2.5PN accelerations of Eq. (|2.13|) . we must include 
the 1PN spin-spin terms in the equations of motion that are generated by the many time derivatives; the explicit 1PN 
corrections to those moments do not contain spin-spin terms. 

The result for the 3.5PN acceleration of body 1 is an expression too lengthy to reproduce here. Calculating a\ — a\ 
and converting all variables to relative coordinates using xi = (m2/m)x and x 2 = — (mi/m)x, we obtain Eq. (| 1 . 3|) . 



B. Spin evolution 

We now want to calculate the evolution of the proper spin S\ to 3.5PN order. A time derivative of Eq. (|2.15[) gives 



S\ = S{[1 



S\ 2vi ■ ai - 3 



m 2 r 



iaKvx-SO-i^ai-Si) 



(4) (4) (3) (3) 

-S\ T' J +S{ X'' -Si 1" +S{ T j 



(3.5) 



We repeat the method of Sec. Ill Bl to determine the contributions of 3.5PN fluid terms to the time derivative of the 
baryonic spin Si, by calculating e yfe J ± p*x 3 a 1 ^ 5PN d 3 x. Only the terms displayed in Eq. (|3.1[) will contribute spin-spin 
terms. Notice that, as we have discussed, the 2.5PN contribution to S\ (actually a spin-orbit term) cancels the first 
two terms in the last line of Eq. (|3.5j) . For ai , which appears in the 1PN terms in Eq. (|3.5j) , we must substitute the 
2.5PN equations of motion; however these make no spin-spin contribution. For S{ in the final 2.5PN terms in Eq. 
(13. 5p we must substitute the 1PN spin-spin precession equations; finally we must use Eq. (|2.15p to convert from Si 
in the 1PN spin-spin terms back to the proper spin S%; the 2.5PN terms there will generate 3.5PN contributions to 
the spin evolution. 

The result is the 1PN spin precession of Eq. (|1.5p . plus a lengthy 3.5PN expression. However, using the fact that, to 
lowest order Si = 0, together with the identities listed in Appendix \E\ it is straightforward to show that our lengthy 
3.5PN expression is almost, but not quite, a pure total time derivative, given by 



(Si) 



3.5PN-SS 



d 

+ ~dl 



(v-S 2 ) + 30f (n-S 2 ) 



(n x Si) 



15H 



(Si x S 2 ) (14 - a)v 2 - (42 - 3a)r 2 + (21 + a) 



in 



-10n[n • (Si x S 2 )] 3(1 + a)v 2 - 15(1 + a)r 2 - (2 - a)- 
-90arn[v • (Si x S 2 )] - 45(1 + 2a)fv[n • (Si x S 2 )] - 5(1 - 8a)v[v • (Si x S 2 )] 
+3 (n • S 2 ) (n x Si) [3(1 + a)v 2 - 15(1 + a)r 2 + (56 + a)— 1 + 27(1 + a)f (v • S 2 ) (n x Si) 



-9(2 - 3a)r (n • S 2 ) (v x Si) + 3(1 - 4a) (v ■ S 2 ) (v x Si, 



(3.6) 



where a = mi/m 2 . The total time derivative can be eliminated by moving it to the left-hand side and absorbing it 
into a redefined proper spin Si, which differs from the original by meaningless 3.5PN correction terms. This is the 
same philosophy by which we absorbed the 2.5PN terms into the initial definition of proper spin in Eq. (|2.15D . In 
the spin-orbit case of Paper III, all the 3.5PN terms could be so absorbed. However, in the spin-spin case, we find 
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that this is not the case, and there is a residual contribution to the spin evolution, given by the first line in Eq. (|3.6[) . 
Because there is no unique way to absorb total time derivatives into Si, the expression for this residual is not unique, 
although its time average over an orbit is. Notice that the residual term is orthogonal to the spin. In other words, 
the magnitude of the proper spin of body 1 is not affected by radiation reaction to 3.5PN order, again not surprising 
for a spinning axisymmetric body. 

Redefining the proper spin, we obtain finally Eq. (|1.6j) . This represents a pure precession of the spin of body 1 
about the radial direction n; however, if the companion spin is perpendicular to the orbital plane, there is no effect. 



C. Comparison with fluxes of energy and angular momentum 

The fluxes of energy and angular momentum in gravitational waves from a binary with spin-spin interactions were 
derived by Kidder et al. [1, [§] , and are given by 



dE 
dt 

~dt 



E N + E ss ■ 
Jn + Jss , 



where we include only the lowest-order "Newtonian" and 1PN spin-spin contributions, given by 

8 r) 2 m 4 



E N 
Ess 

Jjv 



15 r 4 
4 rjrn 2 



(I2v 2 ~ llr 2 ) , 

(Si • S 2 ) (141tj 2 - 165f 2 ) + (n • Si) (n ■ S 2 ) (807f 2 - 504w 2 ) 



5 r 3 
2 ji 



9 S 

' L N ' °- 2 °~ 2 



15 r 6 

71 (v • Si) (v • S 2 ) - 171r (v • Si) (n ■ S 2 ) - 171r (v • S 2 ) (n • Si) 

m 

6 (n • <S 2 ) r - 5 (v • S 2 ) 

r 

6 (n • Si) r — 5 (v • Si 

(771 \ 
18v 2 ~ 30f 2 + 11— j + 6 (v • S 2 j r 

(n • Si) (l%v 2 - 30r 2 + 11—) + 6 (v • Si) r 



ifAi^Si)- 



(n x S 2 ) — 

r 



+ (v x Si) 



(v x S 2 ) 



(Si ■ S 2 ) (60r 2 - 24?j 2 - 46 — ) + 5 (n ■ Si) (n ■ S 2 ) (24tj 2 - 84f 2 + 36- 



-90r ((n • Si) (v ■ S 2 ) + (v • Si) (n ■ S 2 )) - 12 (v Si) (v ■ S 2 



Kidder Q did not determine Jss explicitly, but rather left it in the form of 



(2) 



(3) 



(2) 



(3) 



(2) 



(3) 



(3.7) 
(3.8a) 

(3.8b) 
(3.8c) 



(3.8d) 



Jhs = ~^ 3k { I <il> (a N ) T <kl> (a PN -ss)+ I <jl> (a PN -ss) I <kl> (a N ) + ^ J^ jl> (a N ) J< kl > (a N ) } , (3.9) 



where the mass and current multipole moments are given by Eqs. (|3.2p , where the angular brackets around indices 
denote the symmetric, trace-free part. The notation (aN) or (apN-ss) denotes which acceleration, Newtonian or 
spin-spin, is to be used for the acceleration generated by the time derivatives, and the subscript S denotes the spin 
part of the current moment J . Following this procedure and keeping only terms involving the product of Si with S 2 , 
we obtain Eq. j3~8d)) . 

We now calculate the time derivative of the energy and angular momentum expressions (|2.16l) . and substitute the 
equations of motion, PN spin-spin, 2.5PN point-mass and 3.5PN spin-spin terms, along with the 1PN and 3.5PN spin 
precession equations. After recovering the fact that all 1PN spin-spin contributions cancel, leaving E and J conserved 
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to that order, we find that the changes in E and J due to 2.5PN and 3.5PN spin-spin radiation reaction are obtained 
from the following expressions, 

E = fMV • (a 2 .5PN + a 3 . 5 PN-ss) j 

j = /ix x (a 2 . 5 p N + a 3 . 5PN _ ss ) + S . (3.10) 

Initially, the results do not match the flux expressions above. However, by making use of the identities listed in 
Appendix El we can show that the difference between the expressions in all cases is a total time derivative. These 
can thus be absorbed into meaningless 2.5PN and 3.5PN corrections to the definition of total energy and angular 
momentum. Notice that the residual 3.5PN precession term from S given by the sum of (£1)3. spn and (£2)3. spn horn 
Eq. (|1.6[) exactly balances a corresponding effect in the orbital part, so that the net J matches the flux modulo a total 
time derivative. Thus we have established a proper energy and angular momentum balance between the radiation 
flux and the evolution of the orbit, including spin-spin effects. 



IV. CONCLUSIONS 



We have derived the equations of motion for binary systems of spinning bodies from first principles, including the 
effects of gravitational radiation reaction, and incorporating the contributions of spin-spin coupling at 3.5PN order. 
We found that the spin-spin coupling combined with radiation reaction leads to a small 3.5PN-order precession of 
the individual spins. The resulting equations of motion are instantaneous, dynamical equations, and do not rely on 
assumptions of energy balance, or orbital averaging. They may be used to study the effects of spin-spin interactions 
on the inspiral of compact binaries numerically. We have focusscd attention on effects involving products S1S2 of the 
spins; effects depending quadratically on the individual spins can in principle also be calculated with our approach. 
This will be the subject of future work. 
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APPENDIX A: EXTRACTING TOTAL TIME DERIVATIVES 



Using the Newtonian equations of motion plus the 1PN spin-spin terms, it is straightforward to establish a number 
of identities, which may be used to extract time derivatives from 2.5PN and 3.5PN terms in the expressions (|3 . 6[) and 
(|3.10p . For any non-negative integers s, p and q, we obtain 



d /v 2s f p 



dt V ri 



d /v 2s f p . 



dt V r<i 



v 2s — 2j.p—l 
r q+l 



pv A — (p + q)v 2 f 2 — 2sr 2 



- -pv z - - ip— ((Si ■ S 2 ) - 3 (n • Si) (n • S 2 )) 
r r /i?" 5 



-6s— (f (Si • S2) + (v • Si) (n • S 2 ) + (v • S 2 ) (n • Si) - 5r (n ■ Si) (n • S 2 )) 



v 2s — 2j,p~\ 
r q+l 



4 I , \ 2 -2 o -2 2 111 

pv — (p + q)v r — Zsr pv — 

r r 



3p— ((Si-S 2 )-3(n-Si)(n-S 2 )) 



-6s^r (r (Si • S 2 ) + (v • Si) (n • S 2 ) + (v • S 2 ) (n • Si) - 5r (n • Si) (n • S 2 )) 

/XT" 3 

((n x Si) (n • 5 a ) + (n x S 2 ) (n • Si)) 



(Al) 
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Another set of identities, to be used only in 3.5PN terms, requires only the Newtonian equations of motion: 



dt V ri 



v 2s f p 



dt V r<i 



d /v 2s f p 



dt V r<i 



v 2s-2j,p-l . 

+2v 2 rrx i - l v j) ^ 

v 2s-2f.p-l 
r 9+l I 



— (p + gjw r — 2sf pi) — 

r r . 



2 • 

-2to— ^-x^d- 7 ) 



^23-2^-1 , 



pv 4 — (p + q)v 2 r 2 — 2sr 2 



pv — (p + q)v r — 2sf pv — 

r r . 



pv 



+v rr 



)}■ 



(A2) 



In discussions of inspiralling compact binaries, one often sees the statement that spin-spin terms are 2PN order. This is 
because, for such systems, one treats the radius of each body as being of order m, and the rotational velocity of each 
body V as being of order unity (especially for rapidly rotating black holes); consequently, the spin-spin term in this case 
is effectively of order (m 2 ) 2 /(mr 4 ) ~ (m/r 2 )e 2 . Because the equations derived in this paper apply to arbitrary systems 
treatable with PN methods, we will stick with the formal PN ordering of spin terms. 

A. Papapetrou, Proc. Roy. Soc. (London) 209A, 248 (1951). 

E. Corinaldesi and A. Papapetrou, Proc. Roy. Soc. (London) 209A, 259 (1951). 
G. O'brien, Gen. Relativ. Gravit. 10, 129 (1979). 

B. M. Barker and R. F. O'Connell, Phys. Rev. D 2, 1428 (1970). 
B. M. Barker and R. F. O'Connell, Phys. Rev. D 12, 329 (1975). 

B. M. Barker and R. F. O'Connell, Gen. Relativ. Gravit. 11, 149 (1979). 

L. E. Kidder, C. M. Will and A. G. Wisem an, Phys. Rev. D 47, R4183 (1993) ( |gr-qc/921ld25| . 
L. E. Kidder, Phys. Rev. D 52, 821 (1995) ( |gr-qc/95060"22"j ). 

G. Faye, L. Blanchet and A. Buonanno, Phys. Rev. D 74, 104033 (2006) | gr-qc/0605l39l . 
L. Blanchet, A. Buonanno and G. Faye, Phys . Rev. D 74, 10 4034 (2006) (gr-qc/0605140 ). 

C. M. Will, Phys. Rev. D 71, 084027 (2005) ( |gr-qc/0502039| 

M. E. Pati and C. M. Will, Phys. Rev. D 62, 124015 (2000) (gr-qc/0007087), 
M. E. Pati and C. M. Will, Phys. Rev. D 65, 104008 (2002) (gr-qc/0201001 ). 

See, for example, Sec. 36.8 of C. W. Misner, K. S. Thorne and J. A. Wheeler Gravitation (Freeman, San Francisco, 1973). 
T. Damour and N. Deruelle, Phys. Lett. 87A, 81 (1981). 
T. Damour, C. R. Acad. Sci. Paris 294, 1355 (1982). 



